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Predictive  performance  tools  are  an  important  factor  in  the  success  of  any  new  technology  because  they 
permit  demonstration  of  whether  a  system  will  be  efficient  and  economically  feasible.  A  reliable 
predictive  tool  should  allow  the  designer  to  optimise  system  performance  and  to  maximise  the  cost 
effectiveness  of  a  system  prior  to  installation.  This  article  discusses  the  five  main  parametric  models  for 
photovoltaic  (PV)  systems  available  in  the  literature.  The  current-voltage  characteristics  of  a  photovoltaic 
module  can  be  reproduced  modelling  the  PV  panel  as  an  equivalent  electrical  circuit  made  of  linear  and 
non-linear  components.  The  parameters  describing  such  components  are  directly  related  to  the 
performance  characteristics  of  the  specific  PV  panel,  which  are  generally  available  in  a  graphic  form 
with  respect  to  standard  values  of  temperature  and  incident  irradiance.  Five  of  the  most  recent  models 
are  analysed  in  detail  describing  the  algorithm,  showing  how  the  parameters  were  determined  and 
highlighting  the  most  important  mathematical  and  physical  assumptions  that  characterised  each  model 
and  the  calculation  process.  To  simplify  the  examination  and  reading  of  the  procedures  referred  to  in  the 
five  models,  all  equations  are  rewritten  using  the  same  nomenclature.  Furthermore,  to  assess  the 
accuracy  of  each  model  with  respect  to  the  data  provided  by  the  manufacturer,  this  paper  compares  the 
I-V  and  P-V  curves  at  various  temperatures  and  irradiance  for  a  generic  PV  panel. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


Contents 


1.  Introduction . 685 

2.  Solar  photovoltaic  technologies . 685 

3.  Generic  electric  behaviour  of  a  solid  single  junction  PV  cell . 686 

3.1.  De  Bias,  Torres,  Pieto  and  Garcia  model . 688 

3.1.1.  Mathematical  and  physical  assumptions . 689 

3.1.2.  The  algorithm  steps . 689 

3.2.  Hadj  Arab,  Chenlo  and  Benghanem  model . 689 

3.2.1.  Mathematical  and  physical  assumptions . 690 

3.2.2.  The  algorithm  steps . 690 

3.3.  De  Soto,  Klein  and  Beckman  model . 690 

3.3.1.  Mathematical  and  physical  assumptions . 691 

3.3.2.  The  algorithm  steps . 691 

3.4.  Villalva,  Gazoli  and  Filho  model . 691 

3.4.1.  Mathematical  and  physical  assumptions . 691 

3.4.2.  The  algorithm  steps . 692 

3.5.  Lo  Brano,  Orioli,  Ciulla  and  Di  Gangi  model . 692 

3.5.1.  Mathematical  and  physical  assumptions . 692 

3.5.2.  The  algorithm  steps . 692 

4.  Discussion . 693 

4.1.  Availability  of  necessary  input  data . 693 


*  Corresponding  author.  Tel.:  +39  091  238  619  18;  fax:  +39  091  484  425. 
E-mail  address:  ina@dream.unipa.it  (G.  Ciulla). 

1364-0321/$ -see  front  matter  ©  2014  Elsevier  Ltd.  All  rights  reserved. 
http://dx.doi.Org/10.1016/j.rser.2014.01.027 


C.  Ciulla  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  32  (2014 )  684-696 


685 


4.2.  Synopsis  of  the  major  hypothesis  and  simplifications . 693 

4.3.  Numerical  performance . 693 

4.4.  Overview  of  the  models,  hypothesis  and  simplifications . 695 

5.  Concluding  remarks . 695 

References . 696 


1.  Introduction 

A  steady  increase  in  hydropower  and  the  rapid  expansion  of 
wind  and  solar  power  has  cemented  the  position  of  renewables  as 
an  indispensable  part  of  the  global  energy  mix;  by  2035,  renew¬ 
ables  will  account  for  almost  one-third  of  the  total  electricity 
output.  Solar  power  is  growing  more  rapidly  than  any  other 
renewable  technology.  Renewables  will  become  the  world’s 
second-largest  source  of  power  generation  by  2015  (roughly  half 
that  of  coal),  and,  by  2035,  they  will  rival  coal  as  the  primary 
source  of  global  electricity  [  1  ]. 

In  this  framework,  photovoltaic  generation  systems  have  the 
opportunity  to  be  suitable  for  use  due  to  their  important  advan¬ 
tage  in  being  able  to  produce  electrical  energy  very  close  to  the 
electric  load.  In  this  way,  transmission  losses  are  avoided,  and  it  is 
also  possible  to  satisfy  the  peaks  in  the  daily  load  diagrams 
because  they  supply  the  maximum  power  in  correspondence  to 
the  maximum  request.  Moreover,  photovoltaic  plants  do  not  emit 
pollutant  emissions  or  vibrate;  due  to  their  modularity,  they  can 
comply  with  the  morphology  of  the  installation  sites  [2], 

Promoting  the  dissemination  of  renewable  energy  sources  is 
essential  for  accurate  assessment  of  power  production  to  estimate 
the  economic  and  environmental  performance  with  a  high  grade 
of  reliability  [3,4], 

In  the  field  of  PV  systems,  it  is  essential  to  have  a  good 
simulation  tool  that  can  forecast  the  physical  behaviour  of  a  PV 
panel  by  varying  the  solar  irradiance  and  the  operating  tempera¬ 
ture.  The  operating  temperature  is  also  a  key  variable  of  the 
forecasting  problem  because  it  directly  affects  the  electrical  out¬ 
puts  (current,  voltage  and  power)  and  the  efficiency  of  the  solar 
device.  One  of  the  thermo-physical  processes  that  most  affects  the 
value  of  the  operating  temperature  is  represented  by  convective 
and  radiative  heat  exchange  between  the  panel  and  the  surround¬ 
ing  environment.  In  [5],  it  was  observed  that  the  convective  heat 
transfer  coefficient  is  strongly  influenced  by  the  wind  speed, 
although  it  is  less  affected  by  the  wind  direction  and  is  almost 
unaffected  by  the  air  temperature.  Concerning  the  radiative  heat 
exchange,  the  operating  temperature  also  strongly  depends  on  the 
glass  transparency  and  the  emissivity  and  absorption  coefficient  of 
the  cell.  Some  authors  have  attempted  to  estimate  the  operating 
regimen  of  a  photovoltaic  panel  connected  to  an  electrical  load  in  a 
simple  and  precise  model,  noting  the  operating  temperature  [6-8], 

Some  existing  software,  such  as  PVFORM  [9,10],  PVWATTS  [11], 
PHANTASM  [12],  TRNSYS  [13]  and  PVsyst  [14],  which  is  quite 
sophisticated  and  is  intended  for  advanced  users,  evaluates  the 
physical  behaviour  of  a  photovoltaic  component  and  requires  a 
large  number  of  parameters,  which  are  difficult  to  find  in  technical 
data  and  which  are  often  published  without  the  correlations 
used.  Consequently,  normal  users  usually  employ  the  databases 
included,  sometimes  without  understanding  the  importance  of 
these  data  in  the  assessment  of  the  power  output.  For  these 
reasons,  it  is  certainly  wise  to  learn  more  about  the  features  and 
the  mathematical  approximations  of  the  various  models  used  to 
describe  the  thermo-electrical  behaviour  of  PV  systems. 

The  algorithms  and  models  that  can  be  used  to  assess  the  PV 
power  output  are  generally  of  two  types:  those  that  use  selected 
points  of  the  I-V  curve  [15-19]  and  those  that  exploit  the  entire 


characteristic  curve  [20-23].  The  first  group  of  algorithms  involves 
the  solution  of  five  equations,  using  thermal  drift  coefficients, 
measuring  the  slope  of  the  I-V  curve  and  exploiting  selected 
points,  i.e.,  the  open-circuit,  short-circuit  and  maximum  power 
points.  Although  calculating  the  exact  solution  of  these  equations 
requires  iterative  techniques,  this  method  is  often  much  faster  and 
simpler  compared  to  the  second  group  of  models.  Indeed,  the 
group  of  algorithms  based  on  curve  fitting,  using  deterministic 
optimisation  algorithms,  has  a  disadvantage  that  leads  to  solutions 
for  several  sets  of  local  minima  [24], 

Within  the  first  group  of  algorithms,  the  analytical  models  that 
are  employed  the  most  evaluate  a  solar  cell  as  an  equivalent 
electrical  circuit  composed  of  a  set  of  generally  non-linear  com¬ 
ponents.  The  electrical  parameters  of  the  equivalent  circuit  are 
directly  related  to  the  performance  characteristics  of  the  photo¬ 
voltaic  component,  which  are  generally  described  in  graphical 
form  for  standard  values  of  temperature  and  irradiance. 

The  paper  is  organised  as  follows:  Section  2  presents  the  types 
of  PV,  related  topics  and  their  performances;  and  Section  3 
presents  the  equivalent  circuit  for  PV  modules  and  describes  the 
five  most  recent  and  most  cited  mathematical  models.  An  outline 
of  the  parameter  extraction  procedure  is  given  for  each  model. 
Section  4  presents  a  discussion  of  the  results  obtained  along  with 
the  procedures  used. 


2.  Solar  photovoltaic  technologies 

The  photovoltaic  technologies  are  generally  divided  into  five 

categories: 

1.  Multi-junction  solar  cells :  Multi-junction  solar  cells,  or  tan¬ 
dem  cells,  are  solar  cells  containing  several  p-n  junctions.  Each 
junction  is  tuned  to  a  different  wavelength  of  light,  which 
reduces  one  of  the  largest  inherent  sources  of  losses,  thereby 
increasing  the  efficiency.  Traditional  single-junction  cells  have 
a  maximum  theoretical  efficiency  of  34%;  a  theoretical  “infi¬ 
nite-junction"  cell  would  improve  this  efficiency  to  87%  under 
highly  concentrated  sunlight. 

2.  Single-junction  GaAs :  Gallium  arsenide  is  an  important  optical 
compound  semiconductor  and  is  a  nearly  ideal  material  for 
high  efficiency  single-junction  solar  cells.  In  addition,  GaAs  is 
capable  of  withstanding  radiation  damage  for  a  long  time. 
Owing  to  its  high  efficiency  and  high  resistance  to  radiation, 
GaAs  cells  have  been  favourites  for  powering  satellites  and 
other  spacecraft  [25], 

3.  Crystalline  Si  cells:  Silicon  technology  has  been  the  dominant 
one  for  supplying  power  modules  in  photovoltaic  applications, 
and  the  changes  likely  are  an  increasing  proportion  of  multi¬ 
crystalline  silicon  and  mono-crystalline  silicon  being  used  for 
high-efficiency  solar  cells,  while  the  use  of  thinner  wafers  and 
ribbon  silicon  technology  continues  to  grow  [26,27], 

4.  Thin-film  technologies:  Thin-film  solar  cells  are  basically  thin 
layers  of  semiconductor  materials  applied  to  a  solid  backing 
material.  Thin  films  greatly  reduce  the  amount  of  semiconduc¬ 
tor  material  required  for  each  cell  compared  to  silicon  wafers, 
and  hence,  lower  the  cost  of  production  of  photovoltaic  cells. 
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Nomenclature 

RsO 

reciprocal  of  slope  of  the  I-V  characteristic  of  the 

panel  for  V=V0c  [V]  and  1=0  [12] 

a 

ideality  factor 

Rsh 

shunt  resistance  [12] 

D 

diode  diffusion  factor  (approximately  constant) 

Rsh  0 

reciprocal  of  the  slope  of  the  I-V  characteristic  of  the 

G 

irradiance  [W/m2] 

panel  for  V=V0c  [V]  and  1=0  [12] 

Gsrc 

irradiance  at  SRC  [W/m2] 

SRC 

Standard  Reference  Conditions 

I 

current  generated  by  the  panel  [A] 

T 

cell  temperature  of  the  PV  panel  [K] 

Id 

diode  current  [A] 

V 

voltage  generated  by  the  PV  panel  [V] 

h 

photocurrent  [A] 

vmp 

voltage  at  the  maximum  power  point  [V] 

Imp 

current  at  the  maximum  power  point  [A] 

Yoc 

open  circuit  voltage  of  the  panel  [V] 

he 

short  circuit  current  of  the  panel  [Aj 

lo 

reverse  saturation  current  [A] 

Greek  symbols 

I< 

thermal  correction  factor  [12/  C] 

k 

Boltzmann  constant  [1.381  x  10  23  J/I<] 

a 

ratio  between  the  current  irradiance  and  the 

J 

thermal  correction  factor 

irradiance  at  SRC 

1 VI 

air  mass 

£ref 

material  band-gap  energy  [1.12  eV  for  silicon] 

m 

constant 

Esrc 

amplitude  of  energy  gap  [Jj 

n 

diode  quality  factor 

r 

shape  factor 

Ns 

number  of  cells  wired  in  series 

H'l.sc 

thermal  coefficient  of  the  short  circuit  current  [A/°C] 

<J 

electron  charge  [1.602  x  10  19  C] 

ftv.oc 

thermal  coefficient  of  the  open  circuit  voltage  [V/°C] 

Rl 

electrical  load  [12] 

Rs 

series  resistance  [12] 

GaAs,  copper,  cadmium  telluride  (CdTe),  indium  diselenide 
(CuInSe2)  and  titanium  dioxide  (Ti02)  are  materials  that  have 
been  used  the  most  for  thin-film  PV  cells  [27], 

5.  Other  emerging  photovoltaic:  In  this  category,  there  are  different 
types  of  PV  technology,  such  as  dye-sensitised,  organic,  inorganic 
and  quantum  dot  solar  cells. 

A  dye-sensitised  solar  cell  (DSSC,  DSC  or  DYSC)  is  a  low-cost 
solar  cell  belonging  to  the  group  of  thin-fOilm  solar  cells  [28],  It  is 
based  on  a  semiconductor  formed  between  a  photo-sensitised 
anode  and  an  electrolyte  that  constitutes  a  photo-electrochemical 
system. 

Organic  photovoltaic  (OPV  or  excitonic)  cells  may  be  based  on 
organic  dyes,  semiconducting  polymers,  small  semiconductor 
molecules,  or  on  some  combination  of  these  species.  In  all  cases, 
however,  light  absorption  by  the  organic  species  primarily  results 
in  the  production  of  excitons  rather  than  free  electron-hole  pairs 
[29].  These  cells  are  cheaper  to  produce,  lighter  and  more  flexible 
than  traditional  silicon-based  solar  cells,  and  therefore,  have  a  wider 
range  of  applications. 

Inorganic  cells:  In  thin-film  technologies,  there  exists  a  common 
problem  with  conversion  efficiency  due  to  poor  material  quality; 
photo-generated  electrons  and  holes  cannot  travel  very  far  before 
recombination  (short  free-carrier  diffusion  lengths).  If  a  solar  cell 
can  be  made  using  nanoscale  hetero-junctions,  then  each  photo¬ 
generated  carrier  will  have  less  distance  to  travel,  and  the  problem 
of  recombination  can  be  greatly  reduced.  The  main  factor  for  the 
superior  efficiency  of  inorganic  devices  over  organic  devices  lies  in 
the  high  intrinsic  carrier  mobilities  that  exist  in  inorganic  semi¬ 
conductors.  Higher  carrier  mobility  means  that  charges  are  trans¬ 
ported  to  the  electrodes  more  quickly,  which  reduces  current 
losses  via  recombination  [30], 

Quantum  dot  solar  cells  (QDSCs)  are  based  on  the  Gratzel  cell,  or 
dye-sensitised  solar  cell,  architecture  but  employ  low  band  gap 
semiconductor  nanoparticles  that  possess  such  small  crystallite 
sizes  that  they  form  quantum  dots  (such  as  CdS,  CdSe,  Sb2S3  and 
PbS),  instead  of  organic  or  organometallic  dyes,  as  light  absorbers. 
A  quantum  dot  is  a  portion  of  matter  that  has  electronic  properties 
that  are  intermediate  between  those  of  bulk  semiconductors  and 
those  of  discrete  molecules  [31]. 


In  general,  silicon  wafer-based  photovoltaic  cells  that  include 
single  crystal,  multi-crystalline  and  the  hetero-structure  have 
more  than  80%  of  the  market  share  [32],  The  record  conversion 
efficiency  of  a  multi-crystalline  photovoltaic  cell  is  20.4%,  while  a 
single  crystal  silicon  cell  can  reach  an  efficiency  of  26.4%  (Table  1 ). 

To  evaluate  the  performance  and  price  of  different  PV  module 
technologies  above  summarised,  we  can  refer  to  the  studies 
conducted  by  the  International  Energy  Agency  (IEA)  exemplified 
in  Fig.  1  [33],  or  in  more  recent  studies  [27], 

The  record  conversion  efficiency  of  amorphous  silicon  cells  is 
13.4%,  while  the  record  efficiencies  of  CdTe  and  CuInSe2  are  18.7% 
and  20.4%,  respectively.  Microcrystalline  silicon  has  experienced  a 
breakthrough  recently  and  has  reached  an  efficiency  close  to  that 
of  CdTe.  Third-generation  photovoltaic  cells  are  solar  cells  that 
are  potentially  able  to  overcome  the  Shockley-Queisser  limit  of 
31-41%  power  efficiency  for  single  band-gap  solar  cells.  Common 
third-generation  systems  include  multi-layer  cells  made  of  amor¬ 
phous  silicon  or  gallium  arsenide,  whereas  more  theoretical 
developments  include  frequency  conversion,  hot-carrier  effects 
and  other  multiple-carrier  ejection.  Third-generation  photovoltaic 
cell  technology  is  still  in  the  early  development  stage.  The  cost- 
benefit  ratio  of  this  technology  is  still  too  high  to  be  competitive 
with  the  previous  generation  counterpart. 


3.  Generic  electric  behaviour  of  a  solid  single  junction  PV  cell 

When  a  solid  single  junction  PV  cell  is  not  illuminated,  its 
behaviour  is  similar  to  a  semiconductor  junction;  a  simple  diode, 
which  has  a  characteristic  I-V  curve,  is  traditionally  described  by 
the  equation: 

ID  =  I0(<&- 1)  (1) 

where  /0  is  the  reverse  saturation  current  of  the  diode,  qis  the 
electric  charge  of  an  electron  (1.602  x  10“19  C),  /ds  the  Boltzmann 
constant  (1.381  x  10-23  J/K),  T  is  the  junction  temperature,  VD  is 
the  voltage  across  the  diode  and  n  is  the  ideality  factor,  also  known 
as  the  quality  factor  or  sometimes  the  emission  coefficient.  The 
ideality  factor  n  typically  varies  from  1  to  2  (though  it  can  be 
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higher  in  some  cases),  depending  on  the  fabrication  process  and 
semiconductor  material;  in  many  cases,  n  is  assumed  to  be 
approximately  equal  to  1  (thus,  the  notation  n  often  is  omitted). 

The  illumination  of  the  semiconductor  junction  determines  the 
rise  of  a  photo-generated  current  IL  that  determines  a  vertically 
translation  of  the  I-V  curve  of  a  quantity  linked,  almost  exclusively, 
to  the  surface  density  of  the  incident  energy  (Fig.  2). 

In  this  way,  an  ideal  cell  is  represented  as  a  current  generator, 
connected  in  parallel  with  a  diode,  and  its  I-V  characteristic  is 
described  by  Shockley  [35]  with  the  following  equation: 

I  =  IL-ID  =  IL-Io(f$- 1)  (2) 

where  I  is  the  output  current. 


Table  1 

Best  research-cell  efficiencies  of  2013/03  [34]. 


Category 

Types 

Efficiency  [%] 

Multi-junction  cells 

Four  junction  or  more 

37.80 

(2-teminal  monolithic) 

(no  concentrator) 

Three  junction  (concentrator) 

44.00 

Three  junction  (no  concentrator) 

37.70 

Two  junction  (concentrator) 

32.60 

Two  junction  (no  concentrator) 

30.80 

Single  junction  GaAs 

Concentrator 

29.10 

Thin  film  crystal 

28.80 

Single  crystal 

26.40 

Single  crystal 

27.10 

Thin  film  crystal 

20.10 

Crystalline  Si  cells 

Silicon  heterostructures  (HIT) 

24.70 

Multi-crystalline 

20.40 

Thick  Si  film 

16.00 

Cu  (In,  Ge)  Se2 

20.40 

Thin  film  technologies 

CdTe 

18.70 

Nano-micro-poly-Si 

17.20 

Multi-junction  poly-crystalline 

15 

Amorphous  SiH  (stabilized) 

13.40 

Emerging  photovoltaic 

Dye-sensitized  cells 

11.40 

Organic  cells 

11.10 

Inorganic  cells 

11.10 

Organic  tandem  cells 

12.00 

Quantum  dot  cells 

7.00 

Eq.  (2)  represents  only  a  mere  theoretical  description  because  it 
does  not  take  into  account  the  effect  due  to  the  presence  of  the 
electrodes  both  above  and  below  the  semiconductor  layer  that  are 
necessary  for  collecting  the  charges  that  partially  cover  the  intercept¬ 
ing  surface.  In  [36],  Wolf  correctly  observed  that  in  a  PV  cell,  the 
photocurrent  is  generated  not  only  by  a  single  diode  but  also  by  a 
global  effect  due  to  the  presence  of  a  plurality  of  adjacent  and 
uniformly  distributed  elementary  diodes  along  the  surface  that 
separates  the  two  layers  of  the  semiconductor.  Each  elementary  diode 
is  traversed  by  a  current  that  flows  across  the  semiconductor  layers 
along  a  different  line,  characterised  by  different  electric  resistance  and 
voltage  drop. 

The  transverse  component  of  the  current  IL,  parallel  to  the 
surface  of  the  cell,  should  be  different  for  each  elementary  diode; 
in  this  way,  there  is  a  different  characteristic  I-V  for  each 
elementary  diode.  Because  these  diodes  are  considered  in  parallel, 
the  combination  of  them  determines  the  global  I-V  characteristic 
of  the  PV  cell.  The  transverse  electrical  resistance  is  much  higher 
than  the  electrical  resistance  relative  to  the  direct  component  of  IL. 
The  transverse  component  of  the  current  /,  ,  which  only  occurs  in  a 
real  photovoltaic  cell,  gives  rise  to  a  high  energy  dissipation  that 
drastically  decreases  the  conversion  efficiency  of  the  solar  cell. 

Wolf,  in  [36],  proposed  a  simplified  equivalent  circuit,  shown  in 
Fig.  3,  in  which  there  are  only  a  couple  of  diodes,  one  current  generator 
and  two  resistors  that  take  account  of  the  dissipative  effects  above 
described  and  of  the  presence  of  any  constructive  defects. 

The  resolution  of  the  above  equivalent  circuit  permits  to  obtain 
the  following  expression  of  the  I  current  in  an  implicit  form: 


I  =  IL-1 01 


/  V+1R, 

le  "rt 


V+IRs 

Rsh 


(3) 


where  IL  is  proportional  to  the  irradiance  [35];  f01  and  I02  are  the 
diode  saturation  current  of  the  two  diodes,  Rs  and  Rshare  the  series 
and  parallel  resistances,  respectively;  and  rq  and  n2  are  the  ideality 
factors  of  the  two  diodes.  As  it  is  easy  to  see  from  Fig.  3,  it  is  necessary 
to  know  seven  parameters  ( IL,n n2Jo\-  Io2,Rs ,  ks/i)  to  solve  Eq.  (3). 

In  detail,  the  resistors  Rsh  and  Rs  alter  the  slope  of  the  curve 
before  and  after  the  “knee”  of  the  curve,  respectively,  while  the 
ratio  between  f01  and  I02  changes  the  curvature  [37]. 

Despite  not  being  a  mathematically  indeterminate  problem, 
calculation  of  the  seven  parameters  is  made  difficult  by  the 
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Fig.  1.  Current  performance  and  price  of  different  PV  module  technologies  [33]. 
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Fig.  2.  I-V  Characteristic  of  a  diode  and  illuminated  cell. 


I 


Fig.  3.  Schema  of  simplified  two  diodes  equivalent  circuit. 


implicit  form  of  the  equation  and  by  the  presence  of  two 
exponential  terms.  In  the  scientific  literature,  there  are  few  fully 
explained  models  that  permit  coding  of  the  algorithm  to  retrieve 
the  seven  parameters;  furthermore,  these  models  concern  the 
single  cell  and  are  based  on  particular  assumptions  that  limit  their 
application.  Indeed,  some  authors  admit  that  the  initial  conditions 
strongly  affect  the  resolution  of  Eq.  (3)  [38-40].  In  these  models, 
the  procedures  need  to  be  properly  guided  during  the  initial 
estimation  of  the  parameters  to  avoid  inconsistent  results. 

For  these  reasons,  some  authors  have  preferred  to  employ 
different  correlations  not  based  on  an  electric  model  [41-42], 
whereas  others  have  used  a  simplified  model  based  on  a  single 
diode  [15-19],  as  shown  in  Fig.  4.  This  new  equivalent  circuit  is 
described  by  the  following  implicit  equation: 


I  =  k~ 


') 


V+1RS 

Rsh 


(4) 


in  which  there  are  only  five  unknown  parameters  and  only  one 
exponential  term.  According  to  the  traditional  approach,  the 
photocurrent  IL  depends  on  the  irradiance,  /0  is  affected  by  the 
temperature  of  the  cell  and  n,  Rs  and  Rsh  are  constant. 

Some  authors  have  focused  on  studying  the  one-diode  model 
and  have  suggested  some  directions  for  improvement  and/or 
simplification  that  would  allow  determination  of  the  five  para¬ 
meters  (IL,  n,I0,  Rs,  Rsh)  based  on  the  performance  data  of  the 
modules  commonly  provided  by  the  manufacturers.  Here,  we 
review  five  of  the  best  explained  models  available  in  the  literature 
for  the  evaluation  of  the  electrical  behaviour  of  a  single  junction- 
based  PV  module.  The  models  proposed  by  De  Bias  et  al.  [15] 
(par.  3.1),  Hadj  Arab  et  al.  [16]  (par.  3.2),  De  Soto  et  al.  [17]  (par.  3.3), 
Villalva  et  al.  [18]  (par.  3.4)  and  Lo  Brano  et  al.  [19]  (par.  3.5)  will  be 
described  in  detail. 


I 


Fig.  4.  Schema  of  one-diode  simplified  equivalent  circuit. 


3.3.  De  Bias,  Torres,  Pieto  and  Garda  model 

In  [15],  De  Bias  et  al.  address  the  study  of  the  electrical 
behaviour  of  a  PV  panel  determining  the  five  parameters  solving 
an  equations  system  and  applying  some  modifications  and  sim¬ 
plifications  to  the  traditional  model.  De  Bias  et  al.  state  that  their 
model  permits  to  obtain  a  series  of  I-V  curves  taking  into  account 
the  environmental  variability  with  errors  less  than  1%  in  compar¬ 
ison  with  experimental  results. 

Three  of  the  five  equations  of  the  system,  are  determined 
imposing  that  the  following  points  belong  to  the  l-V  curve, 
expressed  by  Eq.  (4): 

1.  Open  circuit  point;  (V  =  Voc,  I  =  0); 

2.  Short  circuit  point  (V  =  0,  I  =  ISc); 

3.  Maximum  power  point] V  =  Vmp,  I  =  Imp). 


In  these  cases,  Eq.  (4)  becomes 
Voc 


0  =  IL-Io  e^-1  - 


Isc  =  h-Io(e's^-\)- 


Rsh 

hcRs 


‘sh 


(5) 

(6) 

(7) 


/  Vmp+lmpKs  \  Vmp+ImpRs 

Imp=h-k[e  si-  -1J - 

The  substitution  in  Eq.  (6)  of  the  value  of  the  photocurrent  derived 
from  Eq.  (5)  yields  the  following: 


,scl  1+fe)  -'■-'•(«*‘-')-'«p-i) 


Vo 


-la  eV-1 


Isc  l+i 


Rs  \  Vn 


Rsh J  Rsh 


i  /  _ot=  ] 

=  I0  (e^-l-e-^  +  l 


Then, 
Io  = 


Isc  l+i 


Rs\  Voc 


vnr  isr«s 


Rsh )  Rsh 

assuming  that 

voc  'sc% 

e  "T  ►e  "T 

De  Bias  et  al.  hypothesise  that 
Rs  \  Vocl 


1  o  = 


Iscl  1+1 


Rsh)  Rsh 


/esr 


(8) 


(9) 


Similarly,  substituting  Eq.  (7)  the  same  value  of  the  photocurrent, 
the  equation  becomes 

Imp  =  lo  (e^—  l)  +  %* -f„ - 1)  -  (10) 

V  /  Ksh  v  '  Ksh 
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Furthermore,  substituting  Eq.  (9)  in  Eq.  (10)  yields  the  following:  Explicating  Rs  results  in  the  following: 


3 


mp  — 


he 


Voc 

Rsh  . 


1  -e 


Voc-Vn. 

Rsh 


ImpRs 

Rsh 


VqC  -  hcRshO 

nT 


(16) 


or  else 


1+id“,sc 


he  I  1 


Rs 

Rsh 


Voc 

Rsh _ 


Vqc  Voc  —  Vmp 

Rsh  Rsh 

Vmp  +  ImpRs  ~  vocl 

1  -  e - trt - 


After  some  mathematical  manipulation,  it  is  possible  to  state  the 
following: 


_ Vmp+ImpRs~Voc _ 

In  ( ^mp -ISc)0+Rs/Rsh)-(Vmp/Rsh)\ 
111 V  IscO+Rs/Rsh)-(Voc/Rsh)  ) 


(ii) 


The  other  two  conditions  that  permit  solving  the  equation  system 
are  determined  by  imposing  the  equality  of  the  derivative  of 
Eq.  (4)  at  short  circuit  and  open  circuit  points,  respectively,  with 
the  reciprocal  of  the  slopes  of  the  characteristic  curve  in  these 
points;  these  values  are  called  Rsh0  and  Rs0,  and  they  can  be 
graphically  evaluated  from  the  experimental  I-V  curves  issued  by 
the  manufacturers. 


4.  The  derivative  at  the  short  circuit  point  is  as  follows: 


dl 

dV 


o 


^exp^+gL 

l+Rs(^exp^+i 


1 

RshO 


(12) 


The  previous  simplifications  that  have  been  adopted  to  determine 
the  unknown  parameters,  e^r^e^r1  and  are  usually 

fully  satisfied. 

3.3.2.  The  algorithm  steps 

With  this  approach,  adopting  an  iterative  process  that  provides 
the  initial  estimate  of  Rs  value,  the  De  Bias  et  al.  model  retrieves 
the  five  unknown  parameters.  In  detail,  the  algorithm  steps  are  as 
follows: 

•  a  starting  Rs  value  is  hypothesised; 

•  n  and  Rsh  values  are  calculated  from  Eqs.  (11)  and  (15); 

•  Rs  value  is  recalculated  with  Eq.  (16); 

•  an  iterative  procedure  verifies  the  convergence  of  Rs  value; 

•  when  Rs  converges,  the  correct  value  of  Rsh  is  recalculated;  and 

•  the  values  of  IL  and  30.are  calculated  using  Eqs.  (5)  and  (9). 

To  describe  the  behaviour  of  a  PV  module  at  different  condi¬ 
tions  of  temperature  and  radiation,  De  Bias  et  al.  suggest  the 
application  of  the  procedure  described  by  the  International  Stan¬ 
dards  IEC891  [43]  that,  starting  from  the  knowledge  of  current  and 
voltage  values  for  the  given  temperature  and  irradiance,  defines 
the  variation  of  the  I-V  characteristics.  To  adapt  the  saturation 
current  of  the  diode  I0  with  a  variable  temperature  conditions, 
1EC891  uses  the  Townsend  equation  [44]: 


5.  The  derivative  at  the  open  circuit  point  is  as  follows: 


dl_ 

dV 


V  =  Voc 


Rsh 


1 +RS 


e~ar  + 


(13) 


3.1.3.  Mathematical  and  physical  assumptions 

Furthermore,  to  simplify  the  equation  system,  De  Bias  hypothe¬ 
sises  that 


3  o  (scfo  1 
nT6  "T  <Rsh 

so  Eq.  (12)  becomes 

d3_  1  /Rsh  1 

dV  v  =  0  1  +(Rs/RSh)  Rsho 

Finally,  it  is  possible  to  state  that 

Rsh  =  RshO  -  Rs 


(14) 


(15) 


Substituting  Eq.  (9)  into  Eq.  (13)  yields  the  following: 


1 

Rs  0 


'sc(Rsh  +Rs)~vOC 

RshnT  |  1 

nr  'Rsh 


Rsh+Rs- 


?>)■ 


-+1 


Rsh 


and  because  the  value  of  Eq.  (15)  can  be  expressed  as  function  of 
the  resistances  Rsh0  and  Rs0,  the  last  equation  can  be  rewritten  as 
follows: 


l0(T)  =  DT3e& )  (17) 

in  which  eref  is  the  material  band-gap  energy  (1.12  eV  for  silicon  ) 
at  Standard  Reference  Conditions  (SRC),  A  =  y  Ns  where  y  is  the 
completion  or  shape  factor  and  Ns  is  the  number  of  cells  wired  in 
series  (for  ideal  cell  A  is  equal  1)  and  D  is  the  diode  diffusion  factor 
(approximately  constant).  Therefore,  on  the  basis  of  the  correct 
evaluation  of30,  Rs  and  Rsh,  the  other  parameters  can  be  derived 
using  the  values  of  Voc  and  3SC  at  new  conditions  (Eqs.  (5)  and  (6)). 

With  these  assumptions,  the  parameters  that  are  variable  with 
temperature  and/or  irradiance  are  30,  he  and  n.  This  last  parameter, 
which  should  be  considered  independent  of  temperature,  can  be 
evaluated  only  if  the  values  of  Imp  and  Vmp  are  known  in  non¬ 
standard  conditions. 


3.2.  Hadj  Arab,  Chenlo  and  Benghanem  model 


In  [16],  Hadj  Arab  et  al.  proposed  a  five  parameter  model  that  is 
detailed  described  and  experimented  in  [45],  In  this  last  work  the 
electrical  behaviour  of  a  PV  panel  by  referring  to  an  equivalent 
electric  circuit  which  includes  the  photocurrent,  the  reverse 
saturation  current  and  the  series  and  in  parallel  resistances. 
Indeed  the  experimental  test  show  how  the  traditional  model 
with  five  parameters  leads  to  more  accurate  assessments  of 
operating  current  and  more  faithfully  describes  the  behaviour  of 
a  photovoltaic  panel  especially  during  the  maximum  irradiance 
In  this  case,  Eq.  (4)  has  been  rewritten  in  a  different  form: 


,  S’ -ns  , 

I  =  IL-I0  e^r-V 


V  +  IRS 

Rsh 


(18) 


Isc(RShO  -Rs+Rs)-VqC  I  1 
nT  ' 


^s/iO  —  Rs~\~Rs 


I  sc  (Rsho  -Rs  +  Rs)-V  < 

nT 


n  hcRshO-Voe  ,  n  n 


1  + 


+  1] 

hcRsh  0  —Voc 


nT 


in  which  the  five  unknown  parameters  of  the  model  proposed  in 
[25]  are  IL,  30,  Rs,  RSh  and  m  when 


_ [V mp  4~  ImpRsO  —  V oc] _ 

Vt[ln(lsc-^-3mp)  -  ln(/sc-£) 


(19) 
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To  determine  the  values  of  these  parameters,  the  three  known  I-V 
pairs  at  SRC  are  imposed  in  Eq.  (18):  open  circuit  point 
(V  =  V0c ,  1  =  0);  short  circuit  point  (V  =  0,  I  =  /sc)  and  maximum 
power  point  ( V  =  Vmp,  I  =  Imp).  Furthermore  the  resistance  defini¬ 
tions  used  in  Eqs.  (12)  and  (13)  are  applied;  where  Rs0  and  Rsh0  are 
the  inverse  of  the  slope  in  the  open  circuit  point  and  in  the  short- 
circuit  point,  respectively;  these  assumptions  lead  to  a  different 
form  of  Eqs.  (12)  and  (13)  in  which  nT  is  substituted  by  mVt. 


3.2A.  Mathematical  and  physical  assumptions 

To  simplify  the  solution  of  this  system  in  the  model  the 
following  considerations  are  applied: 

If 


1  la 


-e™t 


Rsh  mVt 
Eq.  (13)  becomes 


(20) 


i  IOC 

'0  pmv: 
mVt 

i  +Rstke^ 

So 


1 

RsO 


=  — ;  Ra 


mV, 


.enWt  _  1  + 


mVt 


.gmTf;  Rs  =  R,n  —  ~ 


mVt  voc 
Rs  =  Rs0 — 

*0 

Similarly,  if 
mVt 

Eq.  (12)  can  be  rewritten  as  follows: 

1  /Rsh  1  .  1  1 

1  +(Rs/Rsh)  RshO  Rs+Rsh  Rsl 10 

So, 

Rsh  =  RshO  -  Rs 
Furthermore  if 
Rs<Rsh  —*RshO  =  Rsh 
Similarly,  if 

Voc  „ 
ginVF^l 

in  Eq.  (5),  and  IL  stISc,  then 

Knowing  f0.  Eq.  (7)  can  be  rewritten  as  follows: 


/■Cl /  '  >"P  1  ''/ip-'S  \ 

Imp=Isc-J^(  e  -1  )  - 


Q£- f  Vmp+ImpRs 


,  Voc)  Vmp+ImpRs  -V0c  V mp  (  Rs 
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KsJi  /  Ksh  V  Ksh 


Finally,  if 
Rs 


Rs, 


•  <1 


f  oc 
Jo_plW7 
mVtc 


=  Isc-( Isc-V*)  (eVmp+'^RrVoc - VmP+ImPRs 

V  Rsh  /  V  /  Rsh 

if 

Vmp  +  ImpRs 

e  "Vt  >1 

_  VOC 

then  the  term  e  can  be  neglected: 


(21) 


(22) 


(23) 

(24) 

(25) 

(26) 


(27) 


(28) 


using  the  definition  of  Rs  in  Eq.  (21), 


/  ^OC  1  o 

!SC--n  -  I  e 
Ksh 


Vmp  +lmp  RSQ  ~  vOC 

e  e 


or  rather 

Vmp  +  ImpRs  —  Voc 


1/  ,  i  d  Imp  mV ( 

Vmp  +  lmpKsO - ~  VOC 

•Si - 


=  I  SC  — 


' mp  ,  . 

SC  5  lmp-, 
Ksh 


I  v  mp  I 

,sc~R^~lmP 


-ImpmVt 


mVt  ™Vt(lsc-v1%) 


=  In 


/sC-Vf-ftr 


lSC  R< 


Voc 


(29) 


Thus,  it  is  possible  to  obtain  the  definition  of  mVt  as  explained  in 
Eq.  (19). 

To  solve  the  system  formed  by  Eqs.  (6),  (19),  (24),  (29)  and  (26), 
is  necessary  to  know  the  values  of  fsc  and  Voc.  These  parameters 
are  considered  variables  with  irradiance  and  temperature: 

lsc  =  IscSRCfS - F/q (T  —  TSrc)  (30) 

Csrc 


Voc  =  Voc,sRc+mVt  In 


+  Pvx(T~Tsrc) 


(31) 


3.2.2.  The  algorithm  steps 

The  calculation  procedure  can  be  summarised  as  follows: 

•  the  Rsh  value  is  calculated  by  Eq.  (24); 

•  lsc  and  V0c  parameters  are  calculated  by  Eqs.  (30)  and  (31); 

•  mV t  is  obtained  knowing  Rs0  ,  lsc •  V0c ,  current  and  voltage  at 
the  maximum  power  point  and  using  Eq.  (19); 

•  Rs  is  obtained  with  Eq.  (21);  and 

•  10  is  obtained  from  Eq.  (26); 

•  Finally,  IL  is  obtained  from  Eq.  (6). 

3.3.  De  Soto,  Klein  and  Beckman  model 


De  Soto  et  al.  propose  a  model  in  [17]  for  the  evaluation  of  the 
performance  of  a  PV  panel  represented  by  a  one-diode  equivalent 
circuit  in  which  the  current-voltage  relationship  for  fixed  cell 
temperature  and  irradiance  is  expressed  by  the  following  equation: 


I  =  1,-1, 


y  is.  \  V  +  IRS 

e  a  _l - - - 


(32) 


where  a  is  an  ideality  factor  defined  as  follows: 

NsnkT 

a  =  — - 

Q 


(33) 


Even  in  this  model,  the  solution  of  the  five  parameter  problem 
implies  the  application  of  some  conditions  that  refer  to  the  points 
of  open  circuit  (V  =  V0c,  I  =  0),  short  circuit  ( V  =  0,  /  =  fsc)and  the 
current  and  voltage  at  the  maximum  power  point  ( V  =  Vmp, 

I  =  Imp)  evaluated  at  SRC  The  first  three  equations  of  the  system 
are  as  follows: 


h  —  lo 


\  Voc 
)  R 


(34) 


lsc  =  k-lo(e 


-I  - 


VocRs 


(35) 


/  Vmp  +  Imp  Us  \  Vmp  +  lmp  Rs 

Imp=h-Io[e  «  -1 J - Rh  —  (36) 

where  a  represents  the  ideality  factor  and  IL  is  the  value  of  the 
photocurrent  under  SRC. 
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Because  these  relations  are  insufficient  for  the  determination 
of  the  five  parameters  (IL,  n,  I0,  Rs,  Rsh),  further  equations  are 
obtained  by  imposing  the  derivative  of  the  power  curve  at  the 
maximum  point  equal  to  zero  and  giving  as  known  the  coefficient 
value  of  the  thermal  variation  of  Voc: 


dP 

dV 


P  =  Pmax 


_ .  dl 

—  Imp  +  Vmp-^ 


=  0 


(37) 


The  last  conditions  use  the  expression  of  the  thermal  drift  of  the 
open-circuit  voltage: 


Pv0l 


dV 

'df 


-  —  Vn 


Tsrc  —  T 


(38) 


where  TSrc  is  the  temperature  at  SRC. 

To  determine  /, iVoc ,  it  is  necessary  to  know  the  value  of  the  open 
circuit  voltage  V0c,t •  Because  the  parallel  resistance  Rsh  was 
assumed  to  be  independent  of  temperature,  to  determine  the 
value  of  V0c,T  using  Eq.  (32)  and  then  to  solve  Eq.  (37),  it  is 
necessary  to  know  the  dependency  of  the  parameters  a,  I0  and  IL, 
with  temperature. 


3.3.1.  Mathematical  and  physical  assumptions 

To  this  aim,  the  De  Soto  model  takes  into  account  the  following 
assumptions: 

a)  The  value  of  ideality  factor  a  is  deduced  from  the  temperature 
of  the  cell  T.  In  particular,  for  any  operating  condition,  it  is  a 
linear  function  of  the  temperature  of  the  cell: 

Osrc  Tsrc  ^ 


b)  The  variation  of  the  current  with  the  operating  temperature  is 
explained  by  the  following  equation: 

Iq 

Io.src 

in  which  the  ratio  between  the  current  at  any  conditions  and 
the  current  obtained  at  reference  conditions  I0Src  depends  on 
the  ratio  between  the  temperatures  and  on  an  exponential 
term. 

In  detail,  the  reverse  current  saturation  is  evaluated  with  a 
relation  that  is  derived  from  the  Townsend  equation  Eq.  (17). 
In  Millman-Halkias  [46],  it  is  stated  that 

m  qVC0 

JosKreT  (41) 

where  1<  is  a  constant  and  qVG0  =  sSrc  is  the  amplitude  in  joule 
of  energy  gap;  for  silicon,  the  values  of  r/  and  m  are  2  and  1.5, 
respectively.  Substituting  the  value  of  VT  =  kTNs/q  into  Eq.  (17) 
and  Eq.  (41 ),  it  is  possible  to  state  that  I<  =  D,  m  =  3  ,  r\  =  n. 
Furthermore,  for  an  ideal  cell,  the  equation  becomes 

■3  q  eSRC 

fo ,src  =  D  T SRCek  tsbc  (42) 


T 

Tsrc 


3  1  I  eSRC  _ 


(40) 


the  expression  that  allows  the  determination  of  this  variable  is 
as  follows: 

r — m — [^skc+/</  (T-Tsrc)]  (44) 

C'SRCWlsRC  L  J 


h  = 


where  G  represents  the  irradiance  and  M  is  the  modified  air  mass 
(M  is  an  index  derived  as  function  of  air  mass,  zenith  angle  of  the 
photovoltaic  system  and  of  four  constants  typical  of  each  type  of 
panel;  these  parameters  were  experimentally  obtained  from  NIST 
-  Fanney  et  al.  [47]). 

Finally,  to  complete  the  model,  De  Soto  et  al.  have  analysed  the 
dependence  of  the  resistances  Rs  and  Rsh  with  the  operating 
conditions.  The  Rs  influences  the  shape  of  the  I-V  curves  near  the 
maximum  power  point.  In  [17],  De  Soto  et  al.  show  how,  by 
varying  Rs  by  20%  (in  excess  or  in  deficit)  with  respect  to  the 
reference  value,  the  effect  on  the  curve  is  practically  insignificant, 
thus  justifying  the  assumption  that  Rs  can  be  treated  as  a  constant. 
Rsh  determines  the  slope  of  the  1-V  curve  in  the  short  circuit 
condition,  becoming  more  horizontal  as  it  increases.  In  [17],  the 
authors  claim  that  its  value  seems  to  vary  with  the  irradiance  G  for 
all  types  of  cells;  in  particular,  De  Soto  et  al.  refer  to  the  Schroder 
definition  [48],  which  indicates  that,  for  low  values  of  irradiance, 
the  value  of  Rsh  is  inversely  proportional  to  the  short-circuit 
current  and  therefore  to  the  irradiance.  For  a  description  of  this 
phenomenon,  the  De  Soto  et  al.  model  proposes  the  following 
expression: 


Rsh,SRC  C 

3.3.2.  The  algorithm  steps 

Knowing  the  dependence  of  a ,  I0  and  IL  on  temperature 
(Eqs.  (39),  (40)  and  (44)),  and  of  Rs  and  Rsh  on  the  operating 
conditions,  the  solution  of  the  system  of  five  equations  (Eqs.  (34)-(37) 
and  (38))  with  five  unknowns  is  solved  by  a  nonlinear  solver,  the 
Engineering  Equation  Solver  -  EES  [49], 

However,  the  operating  procedure  adopted  to  solve  this  equa¬ 
tions  system  represents  a  weak  point  in  the  model. 

3.4.  Villalva,  Gazoli  and  Filho  model 

In  [18],  Villalva  et  al.  propose  a  model  to  determine  the 
parameters  of  the  nonlinear  I-V,  shown  in  Eq.  (32),  by  adjusting 
the  curve  at  three  remarkable  points:  open  circuit,  maximum 
power  and  short  circuit  points  provided  by  the  manufacturer. 
Given  these  three  points,  which  are  provided  in  all  commercial 
array  datasheets,  the  method  aims  to  find  the  best  I-V  equation 
including  the  effect  of  the  series  and  parallel  resistances,  and 
warranties  that  the  maximum  power  of  the  model  matches  the 
maximum  power  of  the  real  array. 

3.4.1.  Mathematical  and  physical  assumptions 

Two  other  conditions  are  necessary  to  solve  the  system  of 
equations.  In  [18,50],  the  authors  propose  the  following  assumptions: 


Substituting  Eq.  (42)  into  Eq.  (40),  it  is  possible  to  obtain  the 
relationship  that  describes  the  dependence  of  the  current  /0  on 
the  operating  temperature  of  the  system: 

I0  _  DT3en  _  r  T  l3 
Vsrc  D7"3erw  [Tsrc. 


1  \qJSRC_Q  e 


(43) 


•  The  IL  current  of  the  PV  cell  depends  linearly  on  the  irradiance 
and  the  temperature  according  to  the  following  equation: 

h  =  Ul,  src+Pi  (T-TSrc)]^ —  (46) 

tJSRC 

•  The  authors  propose  a  new  expression  of  I0  : 


c)  The  determination  of  the  photocurrent  IL  is  obtained  as  a 
function  of  short-circuit  current  thermal  drift  plx.  In  particular, 


Isc.src  +Pisc(T  —TSrc) 

vOC£RC  +  PVQC  (T  ~  tSRC ) 

e  ^  - 1 


(47) 


692 


C.  Ciulla  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  32  (2014)  684-696 


Usually,  1  <  a  <  1.5  because  a  expresses  the  degree  of  ideality 
of  the  diode,  and  it  is  totally  empirical;  any  initial  value  can  be 
chosen  to  adjust  the  model.  Eq.  (47)  proposes  a  different  approach 
to  express  the  dependence  of  I0  on  temperature  and  the  thermal 
drift  of  voltage. 


3.4.2.  The  algorithm  steps 

In  [18,50],  the  authors  propose  an  algorithm  for  adjusting  the 
resistance  values  based  on  the  fact  that  there  is  only  one  pair  (Rs, 
Rsh }  that  assures  that  the  calculated  maximum  power  is  equal  to 
the  maximum  experimental  power.  The  authors  claim  that  the 
proposed  model  is  further  improved  by  taking  advantage  of  the 
following  iterative  solution  of  Rs  and  Rsh:  Rs  is  slowly  incremented 
starting  from  Rs  =  0;  the  initial  value  of  Rsh  is  given  by  the 
following: 


D  Vmp  V  OC.SRC  —  Vmp 

Ksh,  min  —  r  j  j 

*  SC, SRC  *mp  Imp 


(48) 


Each  iteration  updates  Rs  and  Rsh  towards  the  best  model  solution. 
The  initial  guess  of  ILSrc  is  obtained  by  the  following  relationship; 

ft, src  =  — —Isc,src  (49) 

Ksh 

The  values  of  Rs  and  Rsh  are  initially  unknown,  but  as  the  solution 
of  the  algorithm  is  refined  along  successive  iterations,  the  values  of 
Rs  and  Rsh  tend  towards  the  best  solution,  and  Eq.  (49)  becomes 
valid  and  effectively  determines  the  light-generated  current  IL 
while  taking  into  account  the  influence  of  the  series  and  parallel 
resistances  of  the  array.  The  software  coding  this  algorithm  is 
available  in  [51], 


3.5.  Lo  Brano,  Orioli,  Ciulla  and  Di  Gangi  model 


Lo  Brano  et  al.  propose,  in  [19],  a  model  for  the  evaluation  of 
the  performance  of  a  PV  panel  represented  by  a  one-diode 
equivalent  circuit  in  which  it  is  postulated  that  Rs,  Rsh  and  /0 
change  with  the  irradiance;  this  hypothesis  is  realistic  because  a 
PV  cell  and  an  illuminated  semiconductor  diode  have  different 
electrical  behaviours.  Based  on  these  assumptions,  a  new  five 
parameter  model  is  proposed,  and  it  is  described  by  the  following 
equation: 


I(ac,  T)  =  acIL(T)-I0(ac,  T ) 


(  "dv+l  'U-rsucIl  +  i  % 

I  e  - 1 


gc[^+7fl7'-rSj;c)]+/ Rs 

Rsh 


(50) 


where  aG  =  G/GSRC  represents  the  ratio  between  the  generic 
irradiance  G  and  the  irradiance  at  SRC  GSrc,  and  IL(T)  is  the 
photocurrent  that  can  be  evaluated  by  the  following  relation: 

Il(T)  =  Il,src+I4Isc(T-TSrc )  (51) 


3.5.1.  Mathematical  and  physical  assumptions 

Although  the  authors  suggest  a  linear  relationship  between  the 
values  of  the  series  and  shunt  resistances  and  the  irradiance,  the 
following  positions  are  made: 


Rs(ac)  =  ^ 
Rs/I(«c)  =  Rt f 


(52) 


Another  important  innovation  concerns  the  determination  of 
the  value  of  I0(aG),  which  can  be  represented  by  the  following 


equation: 


I0(aG,  T)  =  ac 


/lL(T)-Yo<pgry 

I  Ksh 


’ OC(ac.Tl  _ 

e  si  -1 


(53) 


The  variability  of  V0c  with  the  temperature  can  be  evaluated  with 
the  following  equation: 

Voc(aG,T)  =  V0c,src(«  g,  T)+fiVoc(T-TSRc)  (54) 


Lo  Brano  et  al.  observed,  in  [19],  that  for  numerous  commercial 
panels,  there  is  a  regular  dependence  of  the  reverse  saturation 
current  with  the  irradiance  and  temperature  expressed  by  the 
following  formula: 


/0(ac;T)  =  exp 


/ac-200\  J0(l;  T) 

V  l-0.2j  /o(0.2;T) 


In  /o(0.2;  T ) 


(55) 


To  better  reproduce  the  1-V  characteristic  at  generic  temperature  T,  a 
thermal  correction  factor  J  is  defined  as  follows: 


J(ac,  T) 


V mp.SRC  ~  Vmp(CXG,  T ) 

Imp(T  —TSrc) 


(56) 


Here  J  is  a  thermal  correction  factor  similar  to  the  fill  factor  of  the 
curve  described  by  the  1EC891  [43],  To  determine  the  value  of  J,  the 
value  of  T  should  be  chosen  based  on  the  maximum  or  the 
minimum  working  temperature  of  the  PV  module  and  the  max¬ 
imum  range  of  temperature  within  which  the  data  from  the 
manufacturer  are  provided  (in  [19],  a  temperature  of  75  C  is 
considered). 


3.5.2.  The  algorithm  steps 

The  procedure  to  determine  the  five  parameters  requires  a 
preliminary  assessment  of  Rs  and  Rsh.  It  is  also  necessary  to  known 
the  values  of  Isc,src,  Voc, SRcJmP,sRc.  Vmp,sRc  and  the  value  of  Voc  at 
standard  temperature  and  at  the  lowest  irradiance  (typically 
200  W/m2);  all  these  values  can  be  graphically  extracted  from 
the  standard  I-V  characteristic  curve  provided  by  the  manufac¬ 
turer.  To  determine  J  it  is  necessary  to  have  the  values  of  Imp  and 
Vmp  at  different  temperatures  with  respect  to  SRC.  To  determine 
the  values  of  the  unknown  parameters,  the  three  known  I-V  pairs 
at  SRC  are  imposed  on  Eq.  (50),  obtaining  the  same  functional 
relationships  as  the  Bias  et  al.  model  (Eqs.  (5)-(7),  (12)  and  (13)). 

The  system  is  solved  numerically  via  a  double  cycle  of  iterations 
(trial  and  error)  to  retrieve  the  values  of  Rs  and  n.  To  initiate  the 
solving  procedure,  the  following  positions  are  made: 

{h,SRC  =  Isc,src 

n  n  (57) 

Ksh  =  KshO 


which  is  obtained  by  simplifying  Eqs.  (13)  and  (5)  and  taking  into 
account  the  following: 


!Rs<Rsh0 

f0  £RCrs 

Ip,sRce  nTsRC  ^  1 
hTsrc  RshO 


(58) 


The  procedure  comprises  the  following  steps: 


1.  first  attempt  values  of  n  and  of  Rs  are  set; 

2.  Il,src  and  Rsh  are  calculated  with  Eq.  (57); 

3-  Vsrc.  h.sRc  and  Rsh  are  derived  in  sequence,  respectively  from 
Eqs.  (7)’  (6),  and  (12); 

4.  n  is  derived  from  Eq.  (5),  compared  with  the  previous  attempt 
value  and  is  changed  accordingly; 

5.  the  trial  and  error  process  in  steps  3  and  4  is  reiterated  until 
convergence  of  n  is  achieved  with  the  desired  accuracy; 

6.  once  the  first  trial  and  error  process  for  n  has  been  completed, 
Rs  is  derived  from  Eq.  (12),  compared  with  the  previous 
attempt  value  and  changed  accordingly; 
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Table  2 

Synopsis  of  mathematical  relationship  employed  for  each  model. 


Five  parameter  models  /— 1/  relationship  for  1st,  2nd  and  3rd  conditions 


4th  condition 


5th  condition 


De  Bias  et  al. 
Hadj  Arab  et  al. 

Lo  Brano  et  al. 

De  Soto  et  al. 
Villalva  et  al. 


( /  =  lOsc >  0)  =  0 
J  1  =  1(0,  Voc)  =  0 
[l  =  I(Imp,Vmp)  =  0 


I 


V+IRS 

Rsh 


I  =  ach—lo 


/  aG\V  +JHT-Tsrc)]  +  1RS 

(e  °oT 


aG  IV  +JI(T  —  T src)] +IRS 


Rsh 


I  =  k~1 0 


dl  I 

am  =  o 


In  ‘SCp.  . 
jOexp  nr  +Ri- 

Tin  ;SCRS 

l+Rsl^exp  nT  +, 


1 

RshO 


dl  I 

dV\V  =  Voc 


r 


ij 


dPI 

avlp  =  p„ 


=  0 


k  =  [k,SRC  +Hisc  (T  -  Tsrc)]^ 


VVqc 


_  dV\ 

_arl/  =  o  " 

4c^rc  +p/sc  (t  -  rSRC) 


®e(^OC^RC  +PV0C  0  -  ^ SRC)/QW)  


7.  the  first  trial  and  error  process  in  steps  3-5  is  again  started  and 
completely  developed; 

8.  Rs  is  again  derived  from  Eq.  (12)  and  compared  with  the  value 
of  the  previous  attempt  and  changed  accordingly;  and 

9.  the  second  trial  and  error  process  for  Rs  is  reiterated  until  the 
desired  accuracy  is  achieved. 

The  system  is  solved  in  a  non-approximate  form  because  no 
simplified  equations  were  used.  The  above  procedure  is  coded  into 
a  short  software  routine  using  simple  languages.  It  does  not  need 
to  be  guided  towards  solutions  starting  from  initial  fitted  values  of 
the  searched  parameters.  The  procedure  is  more  reliable  because  it 
normalises  all  parameters  with  respect  to  Voc  and  fsc.  Further¬ 
more,  in  [52],  the  Lo  Brano  et  al.  model  was  validated  with  an 
experimental  setup. 


4.  Discussion 

In  this  section,  the  models  described  above  are  compared, 
highlighting  their  peculiarities  while  discussing  the  availability 
of  initial  data  and  the  reliability  of  the  results. 

4.1.  Availability  of  necessary  input  data 

Each  model  is  characterised  by  a  system  of  five  equations  with 
five  unknowns  (1L,  n,  10,  Rs  and  Rsh).  All  considered  models  require 
the  following  three  known  points  of  the  I-V  curve:  open  circuit 
point,  short  circuit  point  and  maximum  power  point.  These  points 
are  always  available  in  the  datasheets  provided  by  the  manufac¬ 
turers,  and  therefore,  the  obtainment  of  these  necessary  values  is 
immediate. 

•  The  De  Bias,  Hadj  Arab  and  Lo  Brano  models  require  the 
measurement  of  the  slope  of  the  l-V  curve  at  selected  points: 
open  and  short  circuit.  Normally,  the  technical  datasheets 
provided  by  the  manufacturer  do  not  include  these  data,  and 
a  graphic  valuation  must  be  carried  out.  This  procedure  is  only 
possible  if  the  datasheets  provide  a  graphical  representation  of 
the  I-V  curve;  however,  a  graphical  estimation  is  always 
approximate,  which  produces  a  measurement  error  that  can 
affect  the  results: 

•  The  De  Soto  model  also  needs  a  graphical  estimation  of  the  I-V 
slope  at  the  maximum  power  point  and  also  requires  the  value 
of  Voc,  the  thermal  drift  coefficient.  Although  this  value  is 
generally  available,  relationship  (38)  is  a  finite  numerical 
approximation  of  a  derivative. 

•  On  the  other  hand,  the  Villalva  model  requires  the  thermal  drift 
coefficients  of  the  open  circuit  voltage  and  short  circuit  current. 
These  values  are  generally  available  in  the  official  datasheets  of 


Table  3 

Mathematical  and  physical  assumptions  for  each  model. 


Assumptions 

De  Bias  et  al. 

Eqs.  (8),  (14),  (17) 

Hadj  Arab  et  al. 

Eqs.  (20),  (22),  (24),  (25),  (27),  (28),  (30),  (31) 

Lo  Brano  et  al. 

Eqs.  (52)— (58) 

De  Soto  et  al. 

Eqs.  (39),  (43),  (44) 

Villalva  et  al. 

Eqs.  (46)-(49) 

commercial  PV  modules.  This  model  does  not  require  any 
graphical  estimation. 


4.2.  Synopsis  of  the  major  hypothesis  and  simplifications 

To  facilitate  the  identification  of  the  mathematical  relationships 
that  permit  coding  of  the  models  into  a  software  package, 
Tables  2  and  3  collect,  for  each  model,  the  five  fundamental 
equations  used  and  their  simplifications.  As  summarised  in  the 
first  column  of  Table  2,  all  the  models  employ  three  remarkable 
known  I-V  pairs  at  SRC.  These  values  are  substituted  in  the 
fundamental  equation;  the  second  column  of  Table  2  collects  for 
each  model  the  mathematical  function  that  links  current  and 
voltage,  I  =  1(1,  V).  The  last  two  columns  gather  the  fourth  and 
the  fifth  conditions,  which  permits  solution  of  the  equation 
system: 

•  The  De  Bias  et  al.  and  the  Hadj  Arab  et  al.  models  are 
characterised  by  the  same  system  of  equations;  however,  the 
different  physical  and  mathematical  assumptions  listed  in 
Table  3  lead  to  different  results.  The  fourth  and  fifth  conditions 
concern  the  graphical  valuation  of  I-V  slope  at  short  circuit  and 
open  circuit  points. 

•  The  first  three  equations  employed  by  the  De  Soto  at  al.  and 
Villalva  et  al.  models  are  identical;  nevertheless,  the  models  use 
the  thermal  drift  coefficients  in  different  ways.  The  De  Soto 
et  al.  model  also  uses  the  slope  of  the  P-V  curve  at  the 
maximum  power  point. 

•  Finally,  the  Lo  Brano  et  al.  model  proposes  a  reformulation  of 
the  general  fundamental  equation  that  directly  takes  into 
account  the  effect  of  temperature  and  irradiance.  The  last  two 
conditions  are  identical  to  those  used  by  the  De  Bias  et  al.  and 
the  Hadj  Arab  et  al.  models. 

4.3.  Numerical  performance 

To  assess  the  accuracy  of  each  model  with  respect  to  the  data 
provided  by  the  manufacturer,  the  I-V  and  P-V  curves  at  vari¬ 
able  temperature  and  irradiance  obtained  by  each  model  are 
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Fig.  5.  I-V  curves  of  the  PV  panel  at  constant  temperature  (25  °C)  and  at  constant  irradiance  (1000  W/m2). 


Fig.  7.  P-V  curves  of  a  PV  panel  at  variable  irradiance  and  constant  temperature. 

represented  in  the  following  graphs  (Figs.  6-9).  In  more  detail,  the 
authors  used  a  specimen  PV  panel  with  the  I-V  curves  drawn  in 
Fig.  5  and  the  tabular  data  collected  in  Table  4. 

Due  to  a  lack  of  information  regarding  a  detailed  calculation 
procedure,  it  was  impossible  to  obtain  any  numerical  results  from 
the  De  Soto  et  al.  model. 

The  graphs  above  show  that  all  of  the  models  describe  the  I-V 
and  P-V  curves  well  when  the  irradiance  is  1000  W/m2  and  the 
cell  temperature  is  25  °C.  By  varying  the  irradiance  and/or  the 
temperature,  the  agreement  with  the  curves  provided  by  the 
manufacturer  changes  considerably.  Defining  the  mean  absolute 


Fig.  8.  l-v  curves  of  PV  panel  at  variable  temperature  and  constant  irradiance. 


Fig.  9.  P  -V  curves  of  a  PV  panel  at  variable  temperature  and  constant  irradiance. 
error  as  follows  (Tables  5  and  6): 

l\/lAV  ^  V1  recalculated  —  manufacturer 
Ivl/ic  2.  n 

i  =  1  1  i, manufacturer 

where  N  is  the  number  of  the  samples,  the  following  tables  show 
the  accuracy  of  each  model  with  respect  to  the  data  issued  by  the 
manufacturer: 

•  Generally,  the  model  that  best  describes  the  I-V  and  the  P-V 
curves  provided  by  the  manufacturer  is  the  Lo  Brano  et  al. 
model,  which  is  characterised  by  an  MAE  of  0.59  for  the  power 
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Table  4 

Specimen  PV  panel  characteristics. 


Voc  IV] 

he  [A] 

v„„„  [V] 

Imax  [A] 

nv„c  [V/  C] 

Vise  [A/  C] 

1000  W/m2  29.35 

8.07 

23.6 

7.567 

-1.07x10- 

1  2.22  xlQ-3 

200  W/m2  27.2 

1.614 

22 

1.508 

Table  5 

Mean  absolute  error  at  constant  temperature  (25  °C)  between  power  data  issued  by 
the  manufacturer  and  each  model  results. 

Mean  average  error 

Global  MAE 

200 

W/m2  K 

400 

W/in2  K 

600 

W/m2  K 

800 

W/m2  K 

1000 
W/m2  K 

De  Bias 

1.36 

1.87 

2.02 

1.48 

0.39 

1.42 

Hadj  Arab 

3.19 

4.82 

5.01 

4.12 

0.39 

3.51 

Lo  Brano 

0.41 

0.68 

1.06 

0.40 

0.39 

0.59 

Villalva 

1.18 

2.33 

2.38 

2.21 

1.70 

1.96 

Table  6 

Mean  absolute  error  at  constant  irradiance  (1000  W/m2) 
issued  by  the  manufacturer  and  each  model  results. 

between  power  data 

Mean  average  error 

25  C  50  C 

75  C 

Global  MAE 

Hadj  Arab 

0.39 

2.07 

3.71 

2.05 

Lo  Brano 

0.39 

0.53 

0.82 

0.58 

Villalva 

1.70 

1.16 

1.96 

1.61 

calculated  at  different  values  of  irradiance  and  constant  tem¬ 
perature  and  by  an  MAE  of  0.58  for  constant  irradiance  and 
different  temperatures. 

•  The  De  Bias  et  al.  and  Hadj  Arab  et  al.  models  tend  to  slightly 
underestimate  the  current  and  power  values  before  the  max¬ 
imum  power  point;  however,  after  the  maximum  power  con¬ 
ditions,  both  models  overestimate  the  currents  and  shift  the 
curve  slightly  to  the  right.  The  De  Bias  model  does  not  allow  for 
calculation  of  the  curve  at  variable  temperatures. 

•  At  the  same  time,  the  Villalva  et  al.  model  tends  to  overestimate 
the  current  and  power  values  under  any  conditions  other 
than  SRC. 


4.4.  Overview  of  the  models,  hypothesis  and  simplifications 

All  examined  models  employed  an  iterative  procedure  to  solve 
a  system  of  five  equations.  However,  each  model  is  characterised 
by  specific  assumptions  and  simplifications  that  inevitably  affect 
the  results: 

•  The  De  Bias  et  al.  model  postulates  that  the  quality  factor  n, 
which  should  be  constant,  varies  with  temperature  and  irra¬ 
diance;  in  general,  the  value  of  n  is  known  only  at  SRC  (Eq. 
(11 )).  In  practice,  for  each  condition  of  irradiance  and  tempera¬ 
ture,  the  model  needs  to  be  recalculated,  estimating  again  the 
values  of  Rs0  and  Rsh0.  To  avoid  uncertainty  in  the  application  of 
corrective  methods,  in  [15],  De  Bias  et  al.  validated  their  model 
using  lSc  and  Voc  values  known  under  specific  conditions  of 
temperature  and  irradiance. 

•  The  Hadj  Arab  et  al.  model  is  tested  in  [16].  This  latest  work  is 
the  one  that  best  describes  the  proposed  model  and  that  best 


matches  the  electrical  behaviour  of  a  solar  panel.  Nevertheless, 
the  model  is  characterised  by  several  simplifications  regarding 
the  exponential  terms  that  could  affect  the  results.  Further¬ 
more,  the  model  requires  that  the  shunt  resistance  be  much 
greater  than  the  series  resistance.  This  condition  is  generally 
always  verified  in  a  high  quality  PV  module. 

•  In  De  Soto  et  al.,  the  five  parameters  are  determined  through  a 
nonlinear  solver,  a  calculation  tool  that  may  have  difficulty  in 
solving  the  equation  system.  Some  tests  performed  with  auto¬ 
matic  solvers  have  not  led  to  automatic  convergence.  In 
addition,  in  the  original  article  [17],  there  is  not  enough 
information  about  the  calculation  steps. 

•  Lo  Brano  et  al.  proposed  a  new  five  parameter  model  that  was 
improved  and  corrected.  The  model  is  characterised  by  some 
significant  differences,  such  as  the  variation  of  /0,  Rs  and  Rsh 
with  irradiance  and  temperature.  Application  of  the  model 
allows  for  obtaining  the  best  results  in  terms  of  comparison 
with  the  experimental  1-V  and  P-V  curves  issued  by 
manufactures. 

•  The  aim  of  the  Villalva  at  al.  model  is  to  fit  the  mathematical  I- 
V  equation  to  the  experimentally  remarkable  point  of  the  I-V 
curve  of  the  practical  array.  A  new  empirical  equation  is 
proposed  to  calculate  the  reverse  diode  saturation  current. 
Furthermore,  an  algorithm  coded  into  a  routine  is  available 
on  the  web  for  public  download. 


5.  Concluding  remarks 

The  /-V  characteristics  of  a  photovoltaic  PV  module  can  be 
reproduced  assuming  a  one-diode  equivalent  circuit  made  of 
linear  and  non-linear  components.  In  this  paper,  five  of  the  most 
recent  mathematical  one-diode  models  for  predicting  the  I-V 
curve  of  a  PV  cell,  which  take  into  account  the  variation  of 
irradiance,  electrical  load  and  operating  temperature,  were  exam¬ 
ined  in  detail.  These  models  reduce  the  problem  to  the  determina¬ 
tion  of  only  five  parameters:  the  photocurrent  IL,  the  diode 
saturation  current  /0,  the  ideality  factor  n,  and  the  series  and 
parallel  resistances  Rs  and  Rsh. 

The  procedure  for  calculating  these  five  electrical  parameters 
constitutes  the  heart  of  the  models  examined,  and  this  work 
allows  for  their  evaluation  with  regard  to  the  fundamental 
hypothesis,  the  availability  of  necessary  data  and  the  reliability 
of  the  results.  For  each  model,  the  fundamental  equations  and  the 
most  important  assumptions  used  were  reported,  and  to  simplify 
the  analysis  of  the  procedures,  all  equations  were  rewritten  with 
the  same  nomenclature  where  possible.  Finally,  an  accuracy 
analysis  was  performed  to  assess  the  capability  of  each  model  to 
reproduce  the  /-V  curve  of  a  specimen  PV  panel.  This  work 
permitted  verification  that  all  the  models  considered  require 
numerical  data  that  are  not  always  available  in  the  official 
datasheets.  Some  models  require  data  obtainable  only  by  applying 
semi-graphic  procedures,  thus  leading  to  an  inevitable  uncertainty 
in  the  measured  values.  A  common  feature  of  all  models  evaluated 
in  this  paper  is  that,  for  solution  of  the  system  of  five  equations, 
they  require  the  application  of  numerical  techniques  that  may  not 
converge  if  the  first  guess  values  are  not  sufficiently  close  to  the 
real  solution. 

Concerning  the  reliability  of  the  results,  all  of  the  models  can 
reproduce  the  I-V  curve  at  SRC  with  sufficient  accuracy,  although 
the  Lo  Brano  et  al.  model  is  slightly  more  accurate  than  the  others. 
However,  although  it  is  not  particularly  accurate  in  the  calculation 
of  the  I-V  and  P-V  curves,  the  Villalva  et  al.  model  does  not  require 
any  graphic  measurements,  and  a  script  is  already  available  on  the 
web. 
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